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(54) Method of synthesizing measurement data of free-form surface 

(57) (a) A surface shape of an object 1 is measured 
from different positions and/or different directions and 
plural partial measurement data 2 including common 
part la is obtained, (bj A parametric surface showing 
the surface sttape of each measurement range is 
formed from the data, (c) A plurality of sampling points 
are set at a predetermined pitch in the common part 
and the Gaussian curvature K and/or an average curva- 
ture H are calculated from the principal curvatures 
(K 1t K 2 ) in the points, (d) Three or more characteristic 
points are selected in accordance with the order from 
larger absolute values of K or H in the common part and 
a normal vector of a plane constructed by three points 
among the points is calculated, (e) the common parts 
are made coincide with each other by moving the para- 
metric surfaces so as to coincide the three characteris- 
tic points with each other and to coincide the directions 
of the normal vectors with each other. Consequently, 
the synthesis data can be formed from the partial meas- 
urement data obtained from the free-form surface. 
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Description 

BACKGROUN D OF THE INVENTION 

5 Technical Field of the Invention 

[0001] The present invention relates to a method of synthesizing shape data by patching partial surface shape data 
of a free-form surface. 

10 Description of the Related Art 

[0002] When the surface shape of a wide object which cannot be measured by one operation or a complicated solid 
three-dimensional object which can be measured only from different directions such as a case where shape data of the 
surface of a mock-up of a vehicle is measured, it is necessary to measure the surface parts of the object from different 
is positions or directions by using a three coordinate measuring machine (digitizer) and to synthesize (patch) partial 
measurement data obtained. 

[0003] As methods of synthesizing the partial measurement data, the following methods are conventionally known. 

(1) Plural partial measurement date is displayed on a CRT and the operator moves the data while seeing the image 
20 so as to coincide overlapped parts with each other (manual method). 

(2) Rotation and translation distances of a measuring machine are preliminarily measured and the measurement 
values are inversely transformed (mechanical measurement method) . 

(3) Markers such as white or light spots are made by ink or the like and positioning is performed (marker method). 

(4) In case of a simple shape such as Z = f(x, y) , undetermined multipliers of its analytic surface are fit by an opti- 
25 mization method such as least-squares method (analytic fitting method). 

[0004] These conventional methods of synthesizing the partial measurement data have, however, the following prob- 
lems. 

[0005] According to the manual method (1), the skill of the operator is required and only synthesis with low accuracy 
30 can be performed. According to the mechanical measurement method (2), since it is necessary to mechanically rotate 

or translate the measuring machine, the size of the measuring machine is increased and the machine is complicated. 

Since an influence of a backlash or the like of the machine is incurred, the method is limited by the mechanical accuracy. 

According to the marker method (3), since markers are indispensable, the method cannot be applied to an object to 

which markers cannot be made. An error can also occur by making the markers. The analytic fitting method (4)has the 
35 feature that the accuracy which is higher than the mechanical accuracy and is almost equal to the measurement data 

can be obtained. An object is, however, limited to a simple shape (analytic shape Z = f(x, y) ). 

SUMMARY OF THE INVENTION 

40 [0006] The present invention is made in order to solve the problems. That is, an object of the invention is to provide 
a method of synthesizing measurement data of a free-form surface, which can form synthesis data with accuracy as 
high as that of measurement data from partial measurement data obtained from a free-form surface which is not limited 
to a simple shape without being limited to mechanical accuracy, without using markers, and without requiring the skill. 
[0007] The inventors of the present invention directed their attentions to a fact such that characteristic points (for 

45 example, maximal points) can be used instead of markers by classifying the local shapes of the free-form surface into 
three types (mountain, valley, and peak) by using the "principal curvature" as an evaluation model which does not 
depend on the coordinate system. The invention is based on the novel ideas. 

[0008] According to the invention, there is provided a method of synthesizing measurement data of a free-form sur- 
face, comprising the steps of: (a) measuring a surface shape of an object from different positions and directions and 

so obtaining plural partial measurement data including a common part; (b) forming a parametric surface showing the sur- 
face shape of each measurement range from the plural partial measurement data; (c) setting a plurality of sampling 
points at a predetermined pitch in the common part on each parametric surface and calculating the Gaussian curvature 
K = k n x k 2 and an average curvature H = (k 1 + k 2 )/2 from the principal curvature (k-, ,k 2 ) in the sampling points, (d) 
selecting three or more characteristic points in accordance with the order from larger absolute values of the Gaussian 

55 curvature K and/or the average curvature H in the common part and calculating a normal vector of a plane constructed 
by three points among the characteristic points; and (e) making the common parts coincide with each other by moving 
the parametric surfaces so as to coincide the three characteristic points with each other and the directions of the normal 
vectors with each other. 
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[0009] The Gaussian curvature K is a product k.,k 2 of the principal curvatures and k 2 of a three-dimensional sur- 
face. It is known that (1) when K > 0, the shape is elliptic, (2) when K = 0, the shape is parabolic, and (3) when K < 0, 
the shape is hyperbolic. 

[0010] The invention uses the characteristic of the principal curvature of a free-form surface. That is, according to the 
5 method of the invention, three or more characteristic points are selected in accordance with an order from larger abso- 
lute values of the Gaussian curvature K and/or the average curvature H in a common part and the characteristic points 
are used instead of markers, thereby making the common parts of parametric surfaces obtained from plural partial 
measurement data coincide with each other and synthesizing (patching) the measurement data with high accuracy. 
[0011] According to a preferable embodiment of the invention, four or more characteristic points are selected in 
io accordance with the order from larger absolute values of the Gaussian curvature K and/or the average curvature H and 
normal vectors of a plurality of planes each constructed by arbitrary three points are calculated in the step (d), and the 
common parts are coincided with each other by moving the parametric surfaces so as to minimize the difference among 
a plurality of combinations of the three characteristic points and the directions of the normal vectors. 
[0012] By the method, the difference between the common parts is minimized by using optimizing means such as the 
75 least-squares method using a plurality of planes each constructed by arbitrary three points and normal vectors, thereby 
enabling the accuracy of synthesis (patching) to be further increased. 

[0013] Preferably, from the principal curvature (k-i ,k 2 ),(1) when k-jSO and k 2 ^0 and k-ik^O, (a) it is determined that 
two curvatures increase. (2) When k^O and k 2 ^0 and k-jk^, (b) it is determined that two curvatures decrease. (3) 
When k i • k 2 < o , (c) it is determined that one of the curvatures increases and the other decreases, (a), (b), and (c) 
20 are labeled as mountain, valley, and peak, respectively, and are displayed by different symbols or colors in an image 
and the parametric surfaces are moved so that the labels coincide with each other. 

[0014] By employing the method at the same time, the coincidence of projections and recesses of the whole common 
parts can be confirmed on a CRT, occurrence of an error due to an influence of noises and the like is avoided, and the 
reliability can be increased. 

25 [001 5] The other objects and advantageous features of the invention will be made clear from the following description 
with reference to the attached drawings. 

PPIEF PESQftlPTlPN OF TH E PRAWINQS 
30 [0016] 

Fig. 1 is an explanatory diagram of a parametric surface; 
Fig. 2 is a flowchart showing the method of the invention; 

Fig. 3(A) is a conceptual diagram showing acquisition of partial measurement data according to the method of the 
35 invention; 

Fig. 3(B) is a conceptual diagram showing formation of parametric surfaces according to the method of the inven- 
tion; and 

Fig. 3(C) is a conceptual diagram showing coinciding operation of common parts of the parametric surfaces 
according to the method of the invention. 



DESCRIPTION OF PREFERRED EMBODIMENTS 

[001 7] The principle of the method of the invention will be first described. 

[0018] Although a free-form surface can be expressed by an implicit function expression f(x, y, z) = 0, a parameter 
45 expression shown by [Expression 1] which is more easily used is employed in the invention. A surface expressed by the 
parameter is called a parametric surface. Fig. 1 is an explanatory diagram of a parametric surface. 



[Expression 1] 
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x = x(ii,v) = y(ti,v) 




55 



where, Cartesian coordinates (x, y, z) on the surface are a differentiable function with respect to parameters u 
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and v and [a, b] shows a rectangle onauv plane. 

[0019] A free-form surface X = X(u, v) can be expressed by the parameters u and v. [Expression 2] is a relational 
expression in differential geometry. The relation is disclosed in, for example, "Curves and Surfaces for Computer Aided 
Geometric Design" (Farin G., 1988, A Practical Guide, Academic Press). 



[Expression 2] 
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JB = E{u,v)=x u x* 
F^F(u J v)^x v x v 
G = G(v,v)=x v x v 

L = L(u,v) = -x tt n u = nx w 

N = iV(u,v) = -x v n v = nx w 



n - " = jr [Xt, A X u j 



||X U A X* 



25 [0020] When it is assumed that X = dv/du , a normal curve k at an arbitrary point X (u, v) is expressed as Expression 
3 in accordance with Expression 2. 

«(x,t) .«(*), * + + 
30 ' E + 2FX + GX* 



[0021 ] The principal curvatures are k-| ,k 2 which are obtained by solving the following Expression 4 and the Gaussian 
35 curvature K is defined as k-|K 2 and the average curvature H is defined as (k + k 2 )/2 . 



40 



45 



LN - M 2 
KlK7= EG-F* 

NE - IMF + LG 
* ,+/C2= EG-F* 



[0022] The projections and recesses of a free-form surface can be defined by using the above-mentioned principal 
curvatures as follows. 

[0023] That is. from the principal curvature (k-, , k 2 ), (1) when k^O and k 2 *0 and k^k^O, (a) it is determined that the 
so two curvatures increase. (2) When k^0 and k ? ^0 and k^k^O, (b) it is determined that the two curvatures decrease. 
(3) When k 1 • k 2 < 0 , (c) it is determined that one of the curvatures increases and the other decreases. Conse- 
quently, for example, by labeling (a), (b), and (c) as mountain, valley, and peak, respectively, and displaying them by dif- 
ferent symbols or colors in an image, the projections and recesses of the free-form surface can be simply and 
accurately grasped by the different symbols or colors. 
55 [0024] According to the method of the invention, the characteristic of the principal curvature is used, the projections 
and recesses of the common parts of the measurement data are classified into mountain, valley, and peak as men- 
tioned above, and larger absolute values of the product and/or sum among them are used as markers. 
[0025] Fig. 2 is a flowchart showing the method of the invention. Figs. 3(A), (B), and (C) are conceptual diagrams. As 
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shown in Fig. 2, the method of the invention comprises five steps (a) to (e). The method will be described hereinbelow 
with reference to Fig. 2 and Figs. 3(A), (B), and (C). 

[0026] As shown in Fig. 3(A), (a) the surface shape of an object 1 is measured from different positions and/or direc- 
tions, thereby obtaining plural partial measurement data 2 (four in the diagram) including common parts la. In this case, 
5 the plural partial measurement data 2 is obtained so as to always include the common part 1a. The partial measure- 
ment data 2 is obtained at the same reduced scale. When the common part 1 a is very narrow (for example, a peripheral 
part or the like), it is preferable to finely set a measurement pitch of the part. 

[0027] The measurement is executed by using, for example, a three coordinate measuring machine (digitizer). In this 
case, the measuring machine is not limited. The machine of a contact type using a stylus, a noncontact type using a 

io laser beam or the like, or other types can be also used. The object 1 can be optionally moved, turned, and the like at 
the time of measurement. In the invention, it is unnecessary to measure the movement/rotation amounts in this case. 
[0028] The object 1 has to be a continuous free-form surface. For example, when the outer surface and the inner sur- 
face are separated, they are measured separately. In this case, it is necessary to input a relative distance (for instance, 
the thickness) between the outer and inner surfaces by another means. 

75 [0029] Subsequently, as shown in Fig. 3(B), (b) parametric surfaces 3. 3' expressing surface shapes of measurement 
ranges are calculated and formed from the plural partial measurement data 2. Although the calculation method is not 
limited as long as it is made by a computer so as to follow the definition of Expression 1 , it is desired that the parametric 
surface is as accurate as measurement data. 

[0030] Subsequently, (c) a plurality of sampling points are set at a predetermined pitch in the common part (part cor- 
20 responding to 1 a) on each of the parametric surfaces 3 and 3' and the Gaussian curvature K and the average curvature 
H are calculated from the principal curvatures (k a , k 2 ) at the sampling points. The calculation can be executed from the 
above Expressions 2 to 4. Preferably, the set pitch of the sampling points is equal to or larger than the measurement 
pitch of the object 1. 

[0031] Further, from the obtained principal curvatures (k t , k^, when (1) and k 2 ^0, and k-ik^O, (a) H is deter- 
25 mined that the two curvatures increase. When(2) k-j^O and k 2 =0 and k-jk^^O, (b) it is determined that the two curva- 
tures decrease. When (3) k 1 • k 2 < 0, (c) it is determined that one of the curvatures increases and the other 
decreases. It is preferable to display an image by labeling (a), (b), and (c) as mountain, valley, and peak, respectively, 
and displaying them by different symbols or colors (for example, brown, blue, and white). 

[0032] Subsequently, (d) three or more characteristic points 4 and 4' are selected in accordance with the order from 
30 larger absolute values of the Gaussian curvature K or the average curvature H in the common parts and normal vectors 
5 and 5' of planes each constructed by three ppints are calculated. Since the partial measurement data 2 is obtained 
at the same reduced scale, the characteristic points 4 and 4' are principally coincided with each other and the normal 
vectors 5 and 5 V are also coincided with each other. Since the partial measurement data 2 is separately obtained, how- 
ever, a slight difference may occur. The difference is at a level equivalent to that of the measurement error, so that the 
35 data is processed on assumption that the difference is substantially the same as the measurement error. 

[0033] Subsequently, as shown in Fig. 3(C), (e) the common parts are coincided with each other by relatively moving 
the parametric surfaces 3 and 3* so as to coincide the three characteristic points 4, 4* with each other and the directions 
of the normal vectors 5, 5* with each other. The calculation and the movement of the parametric surfaces (the three 
degrees of freedom in translation and the three degrees of freedom in rotation) can be easily performed by a known 
40 method with respect to the parametric surfaces expressed by parameters. By the movement, the synthesis (patching) 
of the two partial measurement data 2 whose common parts are coincided is completed. In case of three or more partial 
measurement data 2. by repeatedly executing the above-mentioned processes (especially, (e)), all of the partial meas- 
urement data 2 can be synthesized. 

[0034] As mentioned above, according to the invention, three or more characteristic points 4 are selected in accord- 
45 ance with the order of larger absolute values of the Gaussian curvature K or the average curvature H in the common 
parts and are used as markers, thereby coinciding the common parts of the parametric surfaces 3 obtained from the 
plurality of partial measurement data 2 with each other and enabling the measurement data to be synthesized 
(patched). 

[0035] Preferably, in the step of (d), four or more characteristic points 4 are selected in accordance with the order from 
so larger absolute values of the Gaussian curvature K or the average curvature H and the normal vectors 5 of a plurality 
of pianes each constructed by arbitrary three points are calculated and (e) the common parts are made coincided with 
each other by relatively moving the parametric surfaces 3 so as to minimize the difference in the plurality of combina- 
tions of the three characteristic points 4 and the directions of the normal vectors 5. As the optimization, a known method 
such as the least-squares method or the like can be used. 
55 [0036] As mentioned above, by executing the optimization by using not only the three points in accordance with the 
order from larger absolute values of the Gaussian curvature K or the average curvature H but also the plurality of char- 
acteristic points 4, the accuracy of the synthesis (patching) can be further increased. 

[0037] Although the characteristic points 4 are shown by mountains or valleys in Figs. 3(A), (B), and (C), they can be 
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also shown by peaks. It is preferable to confirm that corresponding characteristic points coincide with each other by 
using not only the magnitude of the absolute value of the Gaussian curvature K or the average curvature H but also the 
symbols or colors (labels) of mountain, valley, and peak. 

[0038] It is also preferable to label the projections and recesses of the surface of the common part as (a) mountain, 
(b) valley, and (c) peak, to display them by different symbols or colors in an image, and after that, by also using the 
above-mentioned method, to move the parametric surfaces so that the whole or most parts of the labels coincide, with 
each other. By using the method at the same time, the coincidence of the projections and recesses of the whole com- 
mon parts can be confirmed on the CRT, the occurrence of an enor due to an influence of noises and the like can be 
avoided, and the reliability can be improved. 

[0039] As mentioned above, by using the method of the invention, like the case of measuring the surface shape data 
of the mock-up of a vehicle, even in case of measuring the surface shape of a wide object which cannot be measured 
by one operation or a complicated three-dimensional object which can be measured only from different directions, the 
synthesized data having the accuracy equivalent to that of the measurement data can be formed from the partial meas- 
urement data obtained from a free-form surface which is not limited to a simple shape without mechanically rotating or 
translating the object, without using markers, and without requiring the skill. 

[0040] By combining the method of the invention with a measuring machine or the like, the whole outward appearance 
can be accurately displayed from partial measurement data obtained by dividedly measuring the parts and the quality 
of the appearance of the body of a vehicle, the face of a person, or the like can be easily and accurately discriminated. 
Also in numerical simulation of the body shape of an aircraft or the like or a functional component such as a turbine 
blade, a junction part can be smoothed and analytic accuracy of fluid flow or the like can be increased. 
[0041 ] That is, the method of the invention can be widely industrially used with a three coordinate measuring machine, 
a CAD system, a CAM system, or a simulation system in the wide-ranged molding processing field or the shape meas- 
uring and evaluating field. 

[0042] Although the invention has been described by some preferred embodiments, it will be understood that the 
range of the right included in the invention is not limited by the embodiments. On the contrary, the range of the right of 
the invention includes all of improvements, modifications, and equivalents included in the scope of the appended 
claims. In case of an object having no characteristic points, the conventional marker method or the like can be also used 
in accordance with necessity. 

Claims 

1 . A method of synthesizing measurement data of a free-form surface, comprising the steps of: 

(a) measuring a surface shape of an object from different positions and/or different directions and obtaining 
plural partial measurement data including a common part; 

(b) forming a parametric surface showing the surface shape of each measurement range from said plural par- 
tial measurement data; 

(c) setting a plurality of sampling points at a predetermined pitch in said common part on each of said paramet- 
ric surfaces and calculating the Gaussian curvature K = k 1 x k 2 and an average curvature H = (k 1 + k 2 )/2 
from the principal curvatures (k-|,k 2 ) in said sampling points, 

(d) selecting three or more characteristic points in accordance with the order from larger absolute values of the 
Gaussian curvature K and/or the average curvature H in the common part and calculating a normal vector of 
a plane constructed by three points among the characteristic points; and 

(e) making the common parts coincide with each other by moving the parametric surfaces so as to coincide the 
three characteristic points with each other and the directions of the normal vectors with each other. 

2. A method according to claim 1 , wherein four or more characteristic points are selected in accordance with the order 
from larger absolute values of the Gaussian curvature K and/or the average curvature H and normal vectors of a 
plurality of planes each constructed by arbitrary three points are calculated in the step (d), and 

further comprising the step (e) of making the common parts coincide with each other by moving the para- 
metric surfaces so as to minimize the difference in a plurality of combinations of said three characteristic points and 
said directions of the normal vectors. 

3. A method according to claim 1 or 2, wherein from said principal curvatures fa, k 2 ), 

(1) when k-j^o and k 2 ^0, and 10^2*0, (a) it is determined that the two curvatures increase, 

(2) when k^0 and k 2 ^0 and k^k^O, (b) it is determined that the two curvatures decrease, and 

(3) when k 1 • k 2 < 0 , (c) it is determined that one of the curvatures increases and the other decreases, and 
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(a), (b), and (c) are labeled as mountain, valley, and peak, respectively, and are displayed by different symbols 
or colors in an image, and the parametric surfaces are moved so that the labels are coincided with each other. 
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Fig. 2 
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(54) Method of synthesizing measurement data of free-form surface 



(57) (a) A surface shape of an object 1 is measured 
from different positions and/or different directions and 
plural partial measurement data 2 including common 
part 1a is obtained, (b) A parametric surface showing 
the surface shape of each measurement range is 
formed from the data, (c) A plurality of sampling points 
are set at a predetermined pitch in the common part and 
the Gaussian curvature K and/or an average curvature 
H are calculated from the principal curvatures (k-, ,k 2 ) in 
the points, (d) Three or more characteristic points are 



selected in accordance with the order from larger abso- 
lute values of K or H in the common part and a normal 
vector of a plane constructed by three points among the 
points is calculated, (e) The common parts are made 
coincide with each other by moving the parametric sur- 
faces so as to coincide the three characteristic points 
with each other and to coincide the directions of the nor- 
mal vectors with each other. Consequently, the synthe- 
sis data can be formed from the partial measurement 
data obtained from the free-form surface. 
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